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Abstract 

In 1912, Otto Sackur and Hugo Tetrode independently developed an equation 
for the absolute entropy of a monoatomic ideal gas. The grand achievement in 
the derivation of this equation was the discretization of phase space for massive 
particles, expressed as 5q5p = h, where q and p are conjugate variables and h 
is Planck's constant. Moreover, Sackur and Tetrode succeeded in corroborating 
their equation by ingeniously utilizing thermodynamic data on mercury. In essence, 
they determined h from the equilibrium of liquid and gaseous phases of mercury 
and compared it with the value previously obtained by Planck from black-body 
radiation. In this paper we discuss the arguments leading to the derivation of the 
Sackur-Tetrode equation and describe also the way how they made their equation 
amenable to an experimental test. Doing a numerical computation similar to that 
of Sackur and Tetrode, we demonstrate that, with modern mercury data, Planck's 
constant can be determined with an accuracy of about one percent. 
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1 Introduction 



The formula for the absolute entropy of a monoatomic ideal gas is named after Otto 
Sackur and Hugo Tetrode who independently derived it in 1912 [H El E]- In classical 
thermodynamics the entropy of a monoatomic ideal gas is 

S(E, V, N) = kN (J In | + In Y. + Sq ) , (1) 

where E, V and iV are the kinetic energy, the volume and the number of atoms, respec- 
tively In classical physics the constant so is undetermined. The achievement of Sackur 
and Tetrode was to compute s . At first sight this does not look very exciting, however, 
in order to compute s$ they had to work out the size of "elementary cells or domains" 
in phase space. Only with this knowledge it is possible to count the number of states in 
classical phase space which is a prerequisite for the computation of Boltzmann's absolute 
entropy given by [U [5] 

S = k\nW. (2) 

In this formula, W is the number of possibilities to realize a system compatible with some 
given boundary conditions. Sackur and Tetrode determined the volume of phase space 
cells as h n where h is Planck's constant and n is the number of degrees of freedom. Until 
then, h was associated mainly with photons, i.e. with massless particles. With the work of 
Sackur and Tetrode it became clear that Planck's constant was relevant also for massive 
particles. In this way h became ubiquitous in statistical physics, more than ten years 
before the advent of quantum mechanics. 

This was an amazing result because a priori Planck's constant in the expression hv 
for the energy of a photon has nothing to do with the phase-space volume associated 
with massive particles. This connection was clarified only later by quantum mechanics. 
We want to stress that the elegance of the work of Sackur and Tetrode derives from the 
combination of theoretical considerations and usage of experimental data with which they 
were able to lend credibility to their result. They did so by successfully applying their 
equation to the then available data on mercury, whose vapor is monoatomic and behaves 
in good approximation as an ideal gas, 

Below we list the articles of Sackur and Tetrode and the achievements therein, written 
in the course of the development of their equation. The titles are literal translations from 
the German titles. 

i. O. Sackur, The application of the kinetic theory of gases to chemical problems [0] 
(received September 1911): In this paper Sackur develops the formula for the en- 
tropy S of a monoatomic ideal g function of the size of the elementary cell. 

ii. O. Sackur, The meaning of the elementary quantum of action for gas theory and 
the computation of the chemical constant pQ (no "received date" , must have been 
written in spring 1912): Here Sackur postulates that the size of the elementary cell 
is h n and obtains the absolute entropy S of a monoatomic ideal gas. Using S, he 
computes the vapor pressure over a solid and makes a comparison with data on neon 
and argon. The numerical results, are, however, not completely satisfying. 
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iii. H. Tetrode, The chemical constant and the elementary quantum of action [2] (re- 
ceived March 7, 1912): Tetrode gives an illuminating derivation of S, assuming that 
the size of the elementary cell is (zh) n . He fits the parameter z by using data on 
the vapor pressure of liquid mercury. Due to some numerical mistakes he obtains 
z « 0.070 

iv. H. Tetrode, erratum to The chemical constant and the elementary quantum of ac- 
tion [2] (received July 17, 1912): Tetrode corrects the numerics and obtains now 
z ~ 1. He acknowledges the papers [U E] of Sackur by noting that the formula for 
S has been developed by both of them at the same time. More precisely, he refers 
to a formula for the so-called "chemical constant" pioneered by Nernst [7j, which 
we will define later. 

v. O. Sackur, The universal meaning of the so-called elementary quantum of action [3] 
(received October 19, 1912): He obtains good agreement (±30%) with the data on 
the vapor pressure of mercury and comments on the paper by Tetrode. 

The paper is organized as follows. In section [2] we describe the different approaches of 
Sackur and Tetrode to derive their equation and add some comments. Since historically 
the corroboration of the Sackur-Tetrode equation by using data on the vapor pressure of 
(liquid) mercury was crucial, we give a detailed account of it in section El Moreover, we 
redo the numerics by using modern mercury data in section H] and obtain a reasonably 
good value of Planck's constant. In section [5] our conclusions are presented. A derivation 
of Kirchhoff's equation, which is used in the numerical computation, is found in the 
appendix. 



2 The Sackur-Tetrode equation 
2.1 Tetrode's derivation 

The starting point of Tetrode's reasoning is the entropy formula (T5]) which should, ac- 
cording to Nernst 's heat theorem [7], give the correct value of the entropy without any 
additive constant. Then he considers a system with n degrees of freedom and phase space 
coordinates qi,...,p n , for which he connects W with the number of configurations of 
phase space points. In order to have a finite entropy, it is necessary to discretize phase 
space, which Tetrode does by introducing "elementary domains" of volume 

5q 1 5p 1 - ■ ■Sq n 6p n = a = (zh) n , (3) 

where h is Planck's constant and z is a dimensionless number. Then he argues that, in 
a system of v identical particles, configurations which are related only by exchange of 
particles should not be counted as different. Therefore, denoting by W the number of 
configurations in phase space, the entropy for such a system is 

w 

S = kln—. (4) 

1 Actually, from Tetrode's equations (12) and (13) we would rather deduce z k, 0.02. 
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This is to avoid the Gibbs paradox and obtain 5* as an extensive quantity, though Tetrode 
does not mention Gibbs in this context. Moving on to the monoatomic gas consisting of 
v = N atoms with mass m and spatial volume V, the number of degrees of freedom is 
n = 3N and, for a given maximal energy E of the gas, the volume occupied in phase 
space is computed by 

V(E,V,N) = Jd 3 x 1 Jd 3 p 1 ---Jd 3 x N Jd 3 p N with ±-(p 2 l+ ...+f N ) <E. (5) 
Utilizing the gamma function, this phase space volume is expressed as 



V 2 

According to the arguments above, the entropy is then given by 

S-kln V{E ' V > N) (7) 
In the last step Stirling's formula is used, to wit the approximations 
IniV! ~ iV(lniV — 1) and lnT ( %- + l) ~ ^ fin 3A 



V 2 J 2 V 2 
for large N. This leads to Tetrode's final result 

for the entropy of a monoatomic ideal gas. 

This derivation is of an amazing lucidity. No wonder that 100 years later it is one of 
the standard methods in modern textbooks. The only amendment to Tetrode's derivation 
comes from quantum mechanics which fixes the size of the elementary domain to h n , i.e. 
requires z — 1; the latter result was obtained by Tetrode through a fit to the data of the 
vapor pressure of mercury. 

From equation (Q with z = 1 we infer that the constant sq of equation (JTJ) is given by 

3 , 4:71171 5 . . 

S »=2 hl ^ + 2' (10) 



2.2 Sackur's derivation 

It is much harder to follow Sackur's line of thoughts. Here we sketch the derivation 
of the entropy formula in [3], because there he gives the most detailed account of his 
derivation. In this paper he first derives Planck's law of radiation by considering a system 
of radiators, before he moves on to the ideal monoatomic gas. In both cases Sackur 
defines a time interval r in which the system is monitored and an energy interval Ae for 
the discretization of energy. For the gas the time r is assumed to be so small that during 
this time collisions between atoms can be neglected. Therefore, during the time interval 



4 



of length r, each of the kinetic energies associated with the three directions in space, e x , 
e y , e z , of every atom can be assumed each to lie in a well-defined energy interval of length 
Ae. In other words, Sackur imagines a three-dimensional energy space with x, y and 
z-axis referring to the kinetic energies e x , e y and e z , respectively, and with energy unit 
Ae on every axis. In this way, the energy space is divided into cubes of volume (Ae) 3 and 
the kinetic energy of every particle lies, during the time interval r, in a well-defined cube. 
If the z-th cube is defined by n^Ae < e^ < (n k + 1) Ae (k = x, y, z) with integers n^, the 
energy £j associated with this cube is, for instance, given by 

Si = (n x + n y + n z ) Ae. (11) 

If in the i-th cube the number of atoms is denoted by Ni, such a number will be the 
smaller the smaller r and Ae are. Consequently, Sackur makes the Ansatz 

N i = Nf{e i ){TAe) 3 , (12) 

which corresponds to the lowest order in an expansion in rAe because there is one factor 
rAe for each axis. In this Ansatz, N is the total number of atoms. The number of 
possibilities to distribute a total of N atoms into r energy cubes is given b}@ 

^ = ^ with N = N l + N 2 + ... + N r . (13) 
According to Boltzmann and Planck, the entropy is obtained by 

S = k\nW = kNlnN - kJ^NilnNi (14) 

i 

for large numbers Ni and the most probable distribution is given by the maximum of S 
under the conditions 

J2 ^ = J2 Nf(ei) (rAe) 3 — N, ]T N e t = £ Nf(e t ) (rAe) 3 e, t = E. (15) 

i i i i 

This procedure resembles the derivation of the canonical ensemble. However, the 
difference is that in this case the total energy of the system is divided into energy intervals 
of length Ae whereas Sackur distributes the energies of single particles among cubes 
in energy space. Finding the maximum of S of equation (fT4l) amounts to computing 
the stationary point of the functional — / de/ln/, under the conditions of a fixed total 
number of atoms and a fixed energy, where the function / is defined in the Ansatz (1121) . 
The sought for stationary point is obtained from the maximum of 

$(/,e) = -/ln/ + (a' + l)/-/3e/, (16) 

where the parameters a' and /3 are Lagrange multipliers: 

|? = - In / + a' - fie = f(e) = e a '~P £ = ae~ Pe with a = e a ' . (17) 



2 Note that the usage of this formula implies that Sackur assumes that the atoms are distinguishable 
because for indistinguishable atoms one would simply have W = 1, for the reason that then every 
decomposition of N into a sum Ni + ■ ■ ■ + N r would correspond to a single state. 
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Eventually, Sackur arrives at the Boltzmann distribution 



Ni = Nae~ p£i (U 



for the number of atoms in an energy box. Plugging this expression into formula ( 1141) for 
the entropy and using equation (fl5"l) . the simple expression 

S = -3kN ln(rAe) - kN In a + k/3E (19) 

for the entropy ensues. 

In equation f|T9|) there are three unknowns: rAe, a and B. At this point, referring to 
Sommerfeld [5], Sackur states that the smallest action that can take place in nature is 
given by Planck's constant h. Therefore, he makes the bold assumption that 

rAe = h, (20) 

which he had already made successfully for the derivation of Planck's law of radiation in 
the same paper. The other two parameters are in principle determined by equation (fl~5l) . 
Sackur then argues that, for simplicity, in the two integrals of equation ( 1T51) summation 
can be replaced by integration. For this purpose he makes the following step: 

2 — 

£ k = ^r- (k = x,y,z) de k = — dp k = — dp k , (21) 

where the x k are the average Cartesian components of the distance covered by the atoms 
during the time r. Then Sackur connects the product of the three average distances with 
the volume V of the gas by equating it with the volume per atom: 

xyz = ^. (22) 

It is hard to understand why this equation should hold, but with equations (|2~T|) and (l22l) 
he effectively introduces an integration d 3 xd 3 p in phase spaceJl Moreover, since Sackur 
nowhere introduces the concept of indistinguishable atoms, he needs the factor 1/N in 
equation (|2"21 for avoiding Gibbs paradox, as we will see shortly. So he ends up with 

V 

t de x de y de z = — dp x dp y dp z (23) 
for the integration in equation (|T5|) and obtains 

1 aVm 3 ( 2tt \ 3/2 , ^ 3aVm 3 ( 2vr \ 3/2 

1 = hv-UJ and E = -iir{-^) ■ (24) 

These two equations are easily solved for a and (3. Plugging the solution 

n 3^ / 3A^ \ 3/2 

8 = and a = — 25 

2N V \4nmEj V ; 



3 These manipulations introduce an ambiguity in the integration boundaries: In de^ the integration is 
from zero to infinity, while in dpk Sackur integrates from minus infinity to plus infinity. 
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into equation (1191) . Sackur arrives at his final result 

ff( ™ = **(§l»f + h4 + § h ^ + §). (26) 

Comparing this expression with Tetrode's result ([9]), we see that there is a difference in 
the last term in parentheses; Sackur has 3/2 while while Tetrode has the correct number 
5/2. Thus 

S(z = l)\ Tetrode -S\ Sackur = kN, (27) 

which Sackur observed and commented upon in [3]. It is interesting to note that in his 
previous paper p]| Sackur actually had the correct number. It is kind of amazing that 
Sackur, with his line of reasoning, arrives at nearly the correct result, being off only by 
kN. This difference is indeed important for the comparison of the entropy formula with 
the data from vapor pressure of mercury [21 13]; anticipating equation (l3~8j) . we see that a 
determination of Planck's constant with Sackur's formula would result in a value which 
is too low by a factor of e -1 / 3 ~ 0.72 where e is Euler's number. 

We conclude this section with a comment on equation (|22|) . We know that S is an 
extensive quantity, i.e. S((E, (V, (N) = (S(E, V, N) holds for all ( > 0. If the factor l/N 
had been absent in equation (|22|) . we would have to replace V by NV in equation (126]) ; 
but then S would not be an extensive quantity, as one can easily check. 



2.3 Discussion 

Let us present here, in particular, for comparison with Sackur's treatment, the derivation 
of the entropy of a monoatomic ideal gas by using the canonical partition function Z. 
Since we are dealing with non-interacting particles, Z is given by 

yN 

Z = W.' (28) 

where Z\ is the partition function of a single particle. The factor 1/N\ is present to take 
into account that the particles are indistinguishable. Then the entropy is given by 

S = k (In Z + PE) = kN (hi ^ + 1 + ?lpj , (29) 

where E is the total energy of the iV particles and p = l/(kT). Furthermore, Stirling's 
formula has been used to replace IniV! by N(lnN — 1). If E/N does not depend on 
N, which is the case for the ideal gas, this equation displays the full dependence on N. 
For the monoatomic ideal gas, in the classical approximation, the single-particle partition 
function is given by the integral 



Z x = -!- / d 3 x [ d 3 p exp [ -P— | = with A 
h 6 Jv J \ 2m J A d 



\/27imkT ^ 

being the thermal de Broglie wave length. The integration domain V is the space taken 
by the gas, i.e. the container with volume V. Plugging Z\ into equation fl29l) yields the 
desired entropy 

S(T,V,N) = kN(\n^ + ^j (31) 
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as a function of temperature, volume and particle number. 

We compare Tetrode's and Sackur's result with the entropy formula (l3Tj) by substi- 
tuting 

E = - NkT (32) 

in equations Q and ( 12T)|) J^ We find what we have announced earlier: Tetrode's result 
exactly agrees with equation f[3"TT) . while Sackur's result differs by kN. We can easily 
locate the origin of the difference. Considering the definitions of a and Z\ and taking into 
account equation ffl~5l) . we find that 

a - jk- (33) 

Insertion of this expression into equation (fT9l leads to the entropy (129]) . with the "1" 
within the parentheses being absent. Effectively Sackur replaces In N\ ~ N(ln N — 1) by 
AHniV in his derivation and does, therefore, not fully take into account indistinguishability 
of the atoms. 

The entropy of the monoatomic ideal g function of the pressure p instead of the 

volume V is obtained with the ideal-gas equation by the substitution V = NkT /p. 

As mentioned in the introduction, Sackur and Tetrode tested their equation on mer- 
cury vapor. This element has seven stable isotopes with various nuclear spins s k [5]. 
Therefore, in principle for mercury one has to add the corresponding residual entropy 

7 

S res (Hg) = Nk £ P k (- lnP fc + \n(2s k + 1)) , (34) 
fc=i 

where the P k are the isotopic abundances (Zlfc-Pfc — 1)> to the Sackur-Tetrode formula. 
Of course, in 1912 the mercury isotopes were not known. However, as we will see in the 
next section, in the mercury test only the entropy difference between gaseous and liquid 
phases is relevant. For both phases, however, the same residual entropy is expected and 
thus iSresCHg) of equation (JM]) drops out. 



3 The vapor pressure of mercury and Planck's con- 
stant 

How to subject the absolute entropy of a monoatomic ideal gas to experimental scrutiny? 
Sackur and Tetrode applied the following procedure. Consider the latent heat L(T) of a 
monoatomic substance for the phase transition from the liquid to the gaseous phase. In 
terms of the absolute molar entropies, the latent heat is given by 

L(T) = T (s vapor (T,p(T)) - s liquid (T,p(T))) , (35) 
4 In Tetrode's formula we set z = 1. 
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where p(T) denotes the pressure along the coexistence curve, i.e. the vapor pressure. If 
the vapor behaves in good approximation like a monoatomic ideal gas, then the Sackur- 
Tetrode equation in the form 

/ ^ 5\ 

s yapoT = R\\n— + ~\ (36) 

with the molar gas constant R can be substituted for s VSupor (T,p(T)). For the liquid phase, 
neglecting the p-dependence, the absolute entropy can be expressed as an integral over 
the heat capacity: 

s**m = / o T dT'^p. (37) 

Note that here the integration includes the solid and liquid phases, and the latent heat of 
melting. After insertion of s vapor and si iqui( j into equation (155]) . one obtains an expression 
for the vapor pressure: 

man = .m + h, (*™m^ + g _ r AT , «<n. m 

yK RT h 3 2 Jo RT v ; 

Similar derivations can be found in pUl EH]. Since equation ( 155j) is a direct consequence 
of equation ( 156"]) . it serves as a testing ground for the Sackur-Tetrode equation. For this 
test not only data on the vapor pressure p(T) are needed, but also data on the latent heat 
L(T) and the heat capacity c p (T) in the condensed phase must be available. While for 
p(T) and L(T) it is sufficient have data in a certain temperature interval, one needs to 
know c p (T) as a function of T down to absolute zero. In 1912 the most comprehensive 
set of data was available on mercury. This was utilized by Sackur and Tetrode to test 
their equation. In this test they followed slightly different approaches. Both employed 
the value of Planck's constant h as determined from black-body radiation and inserted it 
into equation ( 155]) . Then Sackur directly computed the vapor pressure of mercury from 
equation ( 155]) and compared his results with the experimental data, whereas Tetrode 
replaced h in equation ( 138)) by zh and carried out a fit of z to the data. 

Now we want to delineate how Sackur and Tetrode actually performed the numerical 
evaluation of equation (15%]) . We follow the exposition of Sackur in [3] because his account 
is sufficiently detailed and easy to follow. On the right-hand side of equation ( 155]) we have 
to discuss the term with L(T) and the integral. In treating the latent heat as a function of 
T, Sackur uses Kirchhoff's equation — see equation (IA10|) in the appendix. Furthermore, 
he assumes that in the temperature interval he considers, which is from 0° C to 360° C, 
the heat capacity in the liquid phase can be regarded to have the constant value Cp iqmd . If 
at a reference temperature T\ the latent heat is L\, then due to Kirchhoff's equation 

L(T) = U + (J R - cjT id ) (T - Ti). (39) 

The integral on the right-handed side of equation (155]) is treated by splitting it into the 
part in the solid phase, the contribution of the phase transition, and the part in the liquid 
phase. Denoting the latent heat of melting by L m and the melting point by T m , this 
integral reads 

T dr £e£Q = f Tm dr c ^° lld(T/) + h^EA + c w ln Z_. r 4 o) 

T' Jo T> T m p r m K ' 
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Again the approximation that the heat capacity of the liquid is temperature-independent 
has been used. Implicitly the additional approximation that the melting temperature 
T m is independent of the pressure has been made. The final form of the vapor pressure, 
prepared for the numerical evaluation, is thus 

U + fd iquid - \B\ Ti 5 rT m c solid (T') 

lnp(T) = T~>rri +-lnT- / dT' p 1 1 



RT 2 Jo RT' 

L m (Tm) c p qUld ln T i ln (2vrm) 3 / 2 fc 5 / 2 | c^ uid ^ 



RT m R T m h 3 R 

This equation corresponds to Sackur's equation on top of p. 82 of [3] and we have written 
the terms in the same order as there. We have refrained, however, from converting the 
natural logarithm to the logarithm to the base of ten, which was used by Sackur. As 
mentioned earlier, Sackur and Tetrode actually determine the chemical constant, defined 
as 

hTTo x ln k = log l 3 • ^ 

from the data and compare this value of C with the value computed with Planck's constant 
obtained from black-body radiation. At that time, the chemical constant was a commonly 
used quantity. It appears not only in the vapor pressure but also in the law of mass action 
of chemical reactions in the gas phase [7]. Note that the conversion of the logarithm 
mentioned above brings about a division by ln 10 ~ 2.3026 in many places in the equations 
in [21 [3]. 

In equation ( HT1) . in the integral over Cp° lld (T)/T both Sackur and Tetrode use a model 
by Nernst [12] for the specific heat of solid mercury. This model is a kind of Einstein 
model [13] but is sums two frequencies, u and 2u. It is interesting to note that the paper 
of Debye concerning the Debye model [H] has a "received date" July 24, 1912, and is 
thus prior to Sackur's paper [3]. Actually, Sackur refers to it in [3], but only in the part 
concerning Planck's law of radiation; in the integration over the solid phase of mercury 
he uses nevertheless Nernst's model. 

We conclude this section by summarizing and commenting on the approximations 
which lead to equation (j4"Ij) . In essence the following approximations have been made: 



i. The vapor is treated SIS db classical ideal eras. 

ii. The molar volume v i of the liquid is neglected compared to the molar volume v g of 
the vapor. 

iii. In the liquid phase the dependence on p of the isobaric heat capacity is negligible 
in the considered temperature interval. 

iv. There are two technical assumptions which facilitate the numerics: The temperature 
dependence of the heat capacity in the liquid phase is neglected and the melting 
temperature T m is pressure independent. 
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From the first assumption it follows that the heat capacity of a monoatomic vapor is 
constant with the value 

5 



c 



vapor 



R, 



p 2 

which is an important ingredient in equation ( 139|) . The thermal equation of state, 

pV = n m RT, 



(43) 



(44) 



where n m the number of moles of the gas, has been used in equation ( 136]) and in the 
derivation of Kirchhoff 's equation — see appendix. The second assumption, which occurs 
only in the derivation of Kirchhoff 's equation, is well justified because the order of mag- 
nitude of the ratio of the molar volumes is v g /vi ~ 10 3 . To discuss the third assumption 
we note that via the Gibbs potential we obtain the relation 



^liquid 



dp 



-Tv a 2 + 



da 
df 



. . 1 dv 

with a = — — — 

v dT 



(45) 



where a is the thermal expansion coefficient. This equation leads to a linear approximat 
of the heat capacity with respect to the pressure: 



ion 



^liquid 



(T,P) 



„liquid 



{T,p )-T \a 2 + 



da 



v(T,po) (p-po). 



(46) 



P=PO 



The pressure po is a reference pressure. It is well known that the p-dependence of c p for 
liquids is suppressed for two reasons. First of all, the product v p ~ 1 JmoP 1 where v is 
the molar volume of the liquid and p ~ 1 bar is rather small. Secondly, also the thermal 
expansion coefficient a is quite small. Thus, the third assumption is very well justified. 
The assumption of a constant heat capacity for the liquid is not so good. For mercury it 
drops between the melting point at -38.84° C and 200° C by 4% [15]. 



4 Our fit of Planck's constant to mercury data 

It is worthwhile to use the thermodynamic data on mercury available at present and 
employ a slight variation of the method of Sackur and Tetrode described in the previous 
section in order to check the accuracy with which Planck's constant can be determined 
in this way. We follow Tetrode's approach in replacing h by zh in equation ( 138]) . In the 
following we will plug in the modern meanvalue of h and determine z from the data. The 
best modern value of h, recommended by CODATA [16], is 

6.62606957(29) x 10~ 34 Js (47) 

In order to account for the slight temperature dependence of the heat capacity of 
liquid mercury we make the ansatz 

liquid (T) =aQ + ^j, + ^2 (4g) 
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and fit the coefficients a®, a\ and a 2 to the input data from the table presented in [12]. In 
this table one can also read off that from the melting point up to a temperature of about 
200° C the heat capacity of gaseous mercury agrees exactly with the ideal-gas value f|43|) . 
Thus we confine ourselves to the temperature interval from —38.84° C to 200° C, in which 
the ansatz ( 1481) should be sufficient. With equations ( 135]) and ( T4"8l) . and taking into account 
Kirchhoff 's equation, we obtain 

L(T) = L + ^R(T- T ) - a (T - T ) - ± fll (t 2 - T 2 ) - l -a 2 (t 3 - T 3 ) , (49) 

while inserting equation (j4*8]) into the entropy formula (13T|) gives 

T 1 
snquid(T) = s + a In — + «i(T - T ) + -a 2 (T 2 - T 2 ) . (50) 

As a reference temperature we take To = 298.15 K, which allows us to use the enthalpy of 
formation and the standard molar entropy from the CODATA Key Values for Thermo- 
dynamics [T7] : 

T = 61.38 ±0.04 kJmol" 1 , s Q = 75.90 ±0.12 J KT 1 mob 1 . (51) 

The value of sq saves us from the non-trivial task of determining the integral in equa- 
tion (1371) with the boundaries T = and T = T . 

The input data for the vapor pressure of mercury we take from the table in [18]. In 
the legend of this table estimated uncertainties of the vapor pressure values are given, 
which we use in the method of least squares in order to fit the parameter z. A further 
input parameter is the atomic weight of mercury, A = 200.95(2) [9]. The mass value for 
mercury is then m = Au where u is the atomic mass unit. For the determination of h 
from mercury data we can safely neglect errors in the physical constants R, k and u. 

With the above input, our best fit value for z is z = 1.003 at xLm = 4-2. Since we 
have at disposal vapor pressure measurements at 75 temperatures [18] in the considered 
interval, but we determine only one parameter, the number of degrees of freedom is 74. 
For such a large number of degrees of freedom the above value of the minimal x 2 tells us 
that the fit is perfect. We take into account the following sources of uncertainties in z: the 
statistical error determined by x 2 ( z ) — Xmin ± 1> the errors in A, L Q and so, and an error 
in c p . We obtain the uncertainties ±0.0002 for the statistical error and ±0.0005 for the 
error in A. These errors are one order of magnitude smaller than the errors originating 
in L and So which are ±0.004 and ±0.005, respectively. We have no information on the 
error in the heat capacity of liquid mercury in [15J. Therefore, we simply vary a by ±1% 
as a generous error estimate [19]; the resulting uncertainty, however, is smaller than the 
statistical error. In summary, our value of z is 

z = 1.003 ± 0.004 (L ) ± 0.005 (s ). (52) 

Of course, the error estimate above is not a sound statistical computation, but we can 
safely argue that, with existing thermodynamic data on the equilibrium of liquid and 
gaseous phases of mercury, Planck's constant can be determined with an accuracy of about 
one percent. Improving the accuracy of Lq and sq might improve the determination of h, 
but due to the approximations pointed out in the previous section, thermodynamic data 
can most probably never compete with quantum physics data for this purpose. 
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5 Conclusions 



Planck's quantum hypothesis in 1900 was a revolutionary step which he justified by re- 
ferring to Boltzmann, because in this way he could count the number of different photon 
states and compute the entropy of a photon gas by using formula (j2J) • The importance of 
the quantum hypothesis became clear only gradually. In the beginning, Planck's constant 
played a role in loosely connected or seemingly unconnected phenomena. The unified 
perspective was achieved only later with quantum mechanics and quantum field theory. 
However, the importance of the quantum hypothesis for atomic and molecular physics, 
including thermodynamic quantities like heat capacities, was suspected quite early, for 
instance, by Sommerfeld [8] who connected Planck's constant with the "actional in pure 
molecular processes." 

In the beginning, apart from black-body radiation, the phenomena to which the quan- 
tum hypothesis could be applied were scarce. In 1905 Einstein could explain the photo- 
electric effect and in 1907 he put forward the "Einstein model" of the heat capacity of 
solids which could account for deviations from the Dulong-Petit law but gave the wrong 
behavior at low temperatures. This flaw was cured by Debye [H] who developed his 
model practically at the same time as Sackur and Tetrode derived their equation. The 
Bohr model of the atom was to follow in 1913. As a side remark, Ernest Rutherford's 
paper on the atomic nucleus appeared in 1911, in the same year when Heike Kamerlingh 
Onnes discovered superconductivity. For an extensive account of the evolution of the "old 
quantum theory" we refer the reader to [20] . 

Just as Planck more than ten years earlier, Sackur and Tetrode referred to Boltzmann 
in the derivation of their equation. One can view the Sackur-Tetrode equation and its 
successful test with thermodynamic data as one of the very first confirmations of Planck's 
quantum hypothesis. This equation was a quite fundamental step towards modern physics 
as it demonstrated the ubiquity of Planck's constant in statistical physics. We stress 
once more that the outstanding feature of the papers of Sackur and Tetrode was the 
combination of theoretical ideas with an ingenious usage of experimental data. 

One may speculate why the work by Sackur and Tetrode is not that well known in 
the physics community as one would expect from its importance in the development of 
quantum theory and statistical physics. One reason is certainly that both died rather 
young. Sackur (1880-1914), who was actually a physical chemist, died in an explosion in 
the laboratory of Fritz Haber, hardly two years after the Sackur-Tetrode equation. On 
the other hand, Tetrode (1895-1931) was a wunderkind who published his first research 
paper, namely the paper on the Sackur-Tetrode equation, at the age of 17. Later on he 
rather lived in seclusion, though he did publish a few papers which were appreciated by 
the community and kept some contact with eminent contemporary physicists before he 
prematurely died of tuberculosis. 

Acknowledgements: The author thanks E.R. Oberaigner for useful discussions and 
P.O. Ludl for a critical reading of the manuscript. 

Here action means the time integral over the Lagrangian. 
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A Kirchhoff 's equation 



For the transition from the liquid (or solid) to the gaseous phase, this equation relates the 
slope of the latent-heat curve to the difference of the heat capacities across the coexistence 
curve. We derive it here because it is not that well known and also because we want to 
give a fairly self-contained account of the physics around the Sackur-Tetrode equation. 
A derivation of Kirchhoff's equation is also found in |TQl l2Tj - 
The starting point is the equation 

L(T) = T(s 2 (T,p(T))- Sl (T,p(T))) (Al) 

for the molar latent heat of the transition from phase 1 to phase 2. In order to simplify 
the notation we define 



As(T) = S2 (T,p(T))- Sl (T,p(T)), 



(A2) 



and analogously Av and Ac p , the molar volume and heat capacity differences, respectively, 
along the coexistence curve. Taking the derivative of equation ( 1A1I) with respect to T, we 
obtain 



dT dT 



+ T 



OAs 



dp 



dp 
dT' 



(A3) 



Though along the coexistence curve the entropy difference is a function of the temperature 
alone because the pressure along this curve is given by p = p{T), the partial derivatives 
in equation f ]A3j) refer to the original dependence of the entropy on temperature and 
pressure; of course, after performing the derivatives, p has to be replaced by p{T). Next 
we perform three substitutions in equation (IA3I) . Firstly, we note that the molar heat 
capacity is given by 

ds 



T — 

8T 



Secondly, we use the Maxwell relation 

ds 
dp 



dv 



(A4) 



(A5) 



Thirdly, we apply the Clausius-Clapeyron equation 

dp L 
dT ~ TAv' 

With these substitutions equation (IA3j) reads 

dL _ L 1 dAv 

dT~T + Cp ~A^~df~ 



(A6) 



(A7) 



So far this equation is general for a phase transition of first order. Now we argue that, in 
the case of the vapor pressure over a liquid (or solid), the third term on the right-hand 



14 



side of equation ( ]A7j) cancels the first term to a very good approximation. To this end 
we consider the thermal expansion coefficient defined by 



a 



1 dv 
v fff 



(A8) 



which, for an ideal gas, is simply 1/T. Then, using vi/v g <C 1, we derive 



1 dAv 



Av dT 



Olg Ol\ 



1 - 



Vl 



"9 „, 



f 



(A9) 



which proves the cancellation announced abovej^l With this step we finally end up with 
Kirchhoff's equation for the latent heat of vaporization: 



dL 

dT « AC - 



(A10) 



Note that for a liquid far below the critical point usually not only vi <C v g holds but also a\ <C a g . 
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